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Abstract — One way of checking on a pregnant woman is inspection using ultrasonic wave Doppler technique, this check is 
done hy putting Doppler ultrasonic waves over the abdomen to capture the heart beat of the fetus with the addition of the gel 
material, there sebaga in pregnant women who are allergic to the gel material. Alternatively we can use the technique in 
frassonik waves. With a digital design II Chebyshev band pass filter with a frequency range of 2 -3 Hz, the heart rate of 
pregnant women with fetal heart beat can be distinguished. With the help of software MatLab and Xilink will be resolved in 
the manufacture of digital filter terser but. The simulation results with the software quite satisfactory. 
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I. Introduction 

1.1 Background Research 

Since the beginning of pregnancy, it is advisable to conduct regular control to the nearest obstetrician. During the inspection 
control is highly dependent on the doctor or clinic in question, but in general here is given the type of examinations 
performed and it is necessary to ensure the health of mother and baby during pregnancy to child birth. One of these checks is 
the examination with Doppler ultrasound, this check is done by putting Doppler ultrasonic waves over the abdomen to 
capture the heart beat of the fetus with the addition of the gel material, there sebaga in pregnant women who are allergic to 
the gel material. 

1.2 Some researchers History: 

[G.S. Dawes, G.H.A. Visser, J.D.S. Goodman, C.W.G. Redman: 1981] ultrasonic sensors are used, compared to 
ultrasound-based technique with electro cardiography A direct and noted that the short-term index calculation error(STV) 
variability(calculated as the beat-to-beat changes) caused by the approach of ultrasoundreacheslOO%. 

[G.W. Lawson, G.S. Dawes, C.W.G. Redman: 1982] uses ultrasonic sensors to do the same, but for other types of fetal 
monitor, download the error log many of the calculations (200%). In both cases the signal is recorded using a fetal monitor 
from generation before hand, which implement peak detection method to identify fetal heart rate. In further experiments, 
using a fetal monitor with the new generation of the auto correlation function. 

[G.W. Lawson, R. Belcher, G.S. Dawes, C.W.G. Redman : 1983] gained a lot of errors is low on TV ratings: -35%. A 
minus sign means that the pen-ultrasonic approach decreased short-term variability in relation to the FECG. 

[J. Wrobel, J. Jezewski, Dawid Roj, Tomasz P., Robert C., Adam : 2010] Estimated instantaneous reliability of the fetal 
heart rate (FHR) is influenced by autocorrelation techniques commonly used in the channel Doppler ultrasound fetal moni- 
tor today. Has noted greatly reduce the short-term variability, which is quite surprising because as has been shown 
previously, determine fetal monitor FHR signal with satisfactory accuracy in relation to electro cardiography direct reference 
to the fetus. The purpose of this study was to identify the source kesala-han caused by a commonly used approach. Results 
may be possible to developed the method for the correction of FHR variability describes the quantitative index for certain 
types of fetal monitor, which relies on the proposed cancellation of the constant error correction component, which has been 
assigned to an average properties of the autocorrelation function. Although the random error component remaining after 
correction is still not very satisfactory considering the instantaneous values, significanting katan Pen-track assessment of the 
reliability of the measurement of fetal heart rate variability in global cases confirmed one hour 

Fetal heart beat sound (relatively weak), then actually using a sensitive microphone, the fetal heart can be detected by 
strengthening the signal generated by the microphone. The problem is no fetal heart beat in addition there is also another's 
heart beat that results in a stronger heart beat intensity. The frequency of fetal heart rate is higher than the heart beat of the 
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mother (fetal heart rate is 120-160beats per minute (about2-3Hz) while the heart of the mother around 80-90 beats per 
minute), then by technical filtering, then the heart beat the mother can be filtered in order to obtain only the fetal heartbeat. 

Therefore, through this study, researchers will propose to design a digital filter (band pass filter) using approaches in 
frasound (2-3 Hz) are not widely used to monitor fetal heart beat rate, using Simulink and XiLink FPGA tool. 

II. Literature Review 

A filter is a circuit that is designed to skip a frequency band ter-course while weakening all signals outside this band. The 
digital filter is an implementation of the algo-Ritma mathematics into the hardware and / or software that operates on an 
input digital signal for the purpose of filtering is achieved. Jarin-gan-network regular filter is active or passive. Tissues of 
passive filters contain only prisoners, inductor, and capacitor alone. Active filters menggunalan transistor or op -amp plus ta- 
hanan, inductors, and capacitors. But ja-rang inductors used in filters are active, because uku-rannya large and expensive and 
could have Kompo-nen-component in a large persistence. With Filter we can filter out the frequency we want, there are 
several types of filters. As we know, based on the frequency that is passed, in general filter terbagi into 4 yak-ni Low Pass 
Filter; High Pass Filter; Band Pass Filter; and Band Stop Filter. HR filters have infinite impulse response filter where the 
output is a function of input current conditions, previous input and output at an earlier time (causal), better known as 
recursive filters that involve pros-ice feedback and feed forward. For this study, researchers used a type of HR digital filter. 
HR filter is one type of digital filters are widely used in applications of Digital Signal Processing (DSP). HR stands for 
Infinite Impulse Response. Why is it called the re-sponge infinite impulse (infinite)? Due to the feedback in the filter, if you 
are me-enter an impulse (ie a signal T followed by many signals '0'), then on his output will continue to oscillate for their 
feedback, although in practice it will be lost on someday. Advantages include HR filter-kan is a need for fewer coefficient for 
frequency response-uensi steep so as to reduce the amount of computation time. In general, the transfer function of the digital 
filter HR (Infinite Impulse Response) Order N [Smith, 1995] are as follows: 


"00 = ^ = 

A (Z) 


b 0 + b t Z 1 + — + b N Z N 
1 + + — + a N Z~ N 


( 1 ) 


The equation is different is expressed as the equation: 

L L 

Y k=Y b nXk-n ~ ^ Wlc-Ti 


71 = 0 


71=1 


( 2 ) 


Input-output equation HR filter is as follows: 


y(n) = a 0 x(n ) + a x x(n - 1) + 

+a M _ 1 x(n - (M - 1)) 
-b-^yin - 1) - b 2 y(n - 2) 
- b N y(n - (n - N)) 


(3) 


Based HR filter design techniques can be categorized in several design techniques. While the filters FIR (Finite Impulse 
Response) basically have two methods of design, namely the window (the window) and pitajamak the transition band 
(multiband with transition bands). HR digital filter design techniques are based on bilinear transformation from analog 
transfer function prototype. Ana-log transfer function is typically one of the type of transfer function: But-terworth, 
Chebyshev Type 1, Type 2 Chebyshev and Elliptic (or Cauer). The difference between these filter types can be explained by 
looking at the lower pass filter (lowpass) analogous to those described below [Antoni -ou, 1993]. 

a. Butterworth transfer function has escaped under maximum flat magnitude response and magnitude response berku- 
rang horizontally with air-tarnbah frequency. 

b. The transfer function passes under Chebyshev Type I had a magnitude responsive to the ri-ak equilibrium on tape 
passes and a magnitude response is reduced horizontally with increasing frequency in the outer tape passes. 
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c. The transfer function passes under Chebyshev Type II has a magnitude response reduced horizontally on tape passes 
with increasing frequency and magnitude response ripple equilibrium pad prevent a tape. 

d. The transfer function passes under Elliptic mempu-nyai ripple magnitude response of equilibrium on the second 
tape, ribbon or ribbon prevent escapes. 

2.1 Advantages of Digital Filters: 1 

1 . The working area is wide (can be frequency for low and high) with a sampling frequency menga-tour. It can even be 
used in areas Very Low Frequency, where the analog filter is difficult to apply in this frequency area, so it can be 
used for applications elektroni-ka biomedical. 

2. Can provide phase response is absolutely linear. 

3. Because the digital filter implementations can use the programmable pro-cessor, the frequency response can be 
changed automatically. 

4. Some of the input signal (several chan-nel) can be filtered by a digital filter without doubling the hardware for each 
channel. 

5. With the development of technology and fab-rikasi VLSI chip, power consumption and smaller digital filter can 
reduce the cost of production. 

6. Advantages include HR filter requires fewer coefficients for steep frequency response so as to reduce the amount of 
computation time. 

7. In the Order HR filter low frequency response that can produce sharp / steep (almost close to ideal). 

8. The digital filter is very stable against the effects of temperature and time. 

9. The digital filter can handle signals with low fre-keunsi accurate. 

2.2 Digital Filter Chebyshev 

Chebyshev filter is a part of the analog or digital type filter that has a sharper roll-off and more ripple in the passband (type I) 
or stopband ripple (type II) than the Butterworth filter. Chebyshev filter has the property to minimize the error between the 
ideal and the actual characteristics of the filter for the filter range, but the ripple in the passband. This type of filter is named 
in honor of a scientist named Pafnuty Chebyshev, due to their mathematical characteristics derived from the Chebyshev 
polynomials. Due to the passband ripple inherent in Chebyshev filter, a filter that has a smooth response in the passband but a 
more irregular response in the stopband are preferred for some applications [Williams, 2006]. To find or search order (order) 
n on Chebyshev, use the following equation 


n = 


log 


10 g min/10 — 1 
10 g max / 10 — 1 



e = Vl0“ max / 10 - 1 


Ctubyitvev type 1 



Figure. 1. Chebyshev tipe I Lowpass Filter [Ifeachor, 202] 


( 4 ) 

( 5 ) 


Page | 27 


International Journal of Engineering Research & Science (IJOER) 


ISSN: [2395-6992] 


[Vol-2, Issue- 3, March- 2016] 


Chebyshev tppe 2 
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Figure. 2. Chebyshev tipe II Lowpass Filter [Ifeachor, 202] 



2.3 F orm of Direct F orm 


The formulas can be cited in the direct form of mathematical equations as below: 

Y(z ) ( a Q + i^z -1 + a 2 z ~ 2 + — I- a k z _k ) 


H(z) = 


X(z) (1 + ^z 1 + b 2 z 2 + --+h k z fe ) 


2.4 Direct Form I 


H(z) = 


XU **- 1 

1 + S)=i bj-z' 1 


Direct approaches to the realization of the HR filter Direct Form I, where the difference equation is evaluated directly. This 
form is practical for small filter, but may not be efficient and not prkatis (numerically stable) to design intricate or complex 
[If each or, 2002]. In general, this takes the form of the element to delay (delay) 2N (for signal input and output signals) for 
the filter of order N. 


Direct Form Ito its realization can directly use equations quoted as below: 


y(n) = A = ) a^xn-i) - ) aj(yn-j) 

Zj i=0 Zj ;=0 (6) 


From this equation, we can construct a diagram as shown in Figure 3 and Figure 4 


■W 



Figure.3.Realizationof the Direct Form I 
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a2 to 


Figure. 4. Schematic diagram of an HR filter 


2.5 Direct Form II 


Alternatively Direct Form II only requires N delay units, where N is the order of the filter, could potentially have as many as 
half a Direct Form I. This structure is obtained by reversing the order of the numerator and denominator of parts Direct Form 
I, because the two systems are actually linear, and commutative apply for the property. Kemuadian, one will see that there are 
two columns of delay (delay) ie z-1 utilizing the center nets, and these can be combined because there were a lot of repeats. 
The disadvantage is that Direct Form II increases the chances of arithmetic (calculation) over-flow filter or a high Q 
resonance [Ruggiero, 1982]. It has been proven that Q increases, the noise around these two forms of topology increases 
indefinitely. This is because, conceptually first signal passes through the filter all poles (which can usually improve gai at the 
resonance frequency) before generating saturation, then pass through a filter that is all zero (which often weakens as many as 
half of all poles bracing). Realizations of Direct Form II can dikutif in mathematical terms: 


Y(z) = H(z) . X(z) 


aw = 


D(z) 


¥ 0) 


NfczJ 

D(z) 


X(z) 


X{z) 

n(z) 


■*(*) 


W{z) = 


xfrO 

D(z) 


Y(z) = N(z) . W(z) 

ivGO = — = — - 
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( 1 F i^z -1 F--+ i^z - *}. M/(z) = X(z) 
w-’(n) -f i L (ivn — 1) -f ■" -h ijtfwra — k) = 
w(n) = x(n) — tqCu-n — 1) — — — — ft) 


ft 


u- (n) = x(n) — 
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(7) 


Y(z) = N(z) . W(z) 

¥(z) = Co E + ^z' 1 + — + a- K z~ K X IVtz] 


;y(n) = %(urn) -f a L ( w n — 1) -f — |- — ft) 
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L = D 


From the equation above, we can construct a diagram as in Figureb below: 


a 0 



Figures. Realization of theDirect Form II 


2.6 Research Objectives 

The specific objectives of this program are: 

1. Creating a digital filter design type of Band Pass Filters HR Chebyshev II order 4 so that it candiimplcment as ikan 
on the FPGA board. 

2. Realizing a digital infrasonic filter that can later be used for development, detection of fetal heart beat in the form of 
SoC (System on Chip). 

3. Program MATLAB is used to memu-dahkan seek coefficient digital filter design. 

4. Program Xlink is used for the simulation of digital filter design MatLab and to mens in kronasikan (match) results. 

2.7 Research Methods 

2.7.1 Digital Filter Design 

This section discusses the design, designing band pass Chebyshev digital filter by performing mathematical calculations. To 
be implemented in a Xilinxlsel4.2 with VHDL is a hardware description language standard for VHSIC. Hardware 
description language used in electronic design automation to describe digital and mixed -signal systems such as the field 
programmable gate array and integrated circuits. Block diagram of the overall system can be described as below: 



Figure 6. Block Diagram Global Research 
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The location is on the pole Butterworth filter and the Chebyshev polynomial coefficients appropriate and have been obtained 
and tabulated by Weinberg (see Table 1). 

An abbreviated set of coefficients polynomial for Butterworth and Chebyshev functions presented here, some form of 
polynomial factors have been calculated and tabulated. Polynomial given in the foundation needed to build a low-pass 
transfer function. The general form of the transfer functions as in the table below: 

TableI. 

The coefficient of low-pass filter polynomial denominator (denominator polynomials) 


Order 

B0 

B1 

B2 

B3 

B4 

B5 




BUTTERWORTH 




1 

1 

1 





2 

1 

1.414 

1 




3 

1 

2 

2 

1 



4 

1 

2.613 

3.414 

2.613 

1 


5 

1 

3.236 

5.236 

5.236 

3.236 

1 



CHE BYSHEV0.5 DB RIPPLE 



1 

2.863 

1 





2 

1.516 

1.425 

1 




3 

0.716 

1.535 

1.253 

1 



4 

0.379 

1.025 

1.717 

1.197 

1 


5 

0.179 

. 

1.309 

1.937 

1.174 

1 



CHEBYSHEV 1 DB RIPPLE 



1 

1.965 

1 





2 

1.103 

1.098 

1 




3 

0.491 

1.238 

0.988 

1 



4 

0.276 

0.743 

1.454 

0.953 

1 


5 

0.123 

0.581 

0.974 

1.689 

0.937 

1 



CHEBYSHEV 2 DB RIPPLE 



1 

1.308 

1 





2 

0.823 

0.804 

1 




3 

0.327 

1.022 

0.738 

1 



4 

0.206 

0.517 

1.256 

0.716 

1 


5 

0.082 

0.459 

0.693 

1.499 

0.706 

1 



CHEBYSHEV 3 DB RIPPLE 



1 

1.002 

1 





2 

0.708 

0.645 

1 




3 

0.251 

0.928 

0.597 

1 



4 

0.177 

0.405 

1.169 

0.582 

1 


5 

0.063 

0.408 

0.549 

1.415 

0 574 

1 


G(s) — 


^0 

( Bq + B^s + B^s 2 ... + Bfcsk') 


( 8 ) 


The coefficient of the denominator of the equation (8) are tabulated in Table to function for the Butterworth and Chebyshev 
functions to four possible values of passband ripple (0.5 dB, 1 dB, 2 dB and 3 dB). The use of the filter order from k = 1 to k 
= 5. If desired gain in dc numerator constant then have AO = BO (Note: the coefficient is defined according to the equation 
( 8 ) ( 6 ). 


Creating a Digital Design Band-Pass Filter Chebyshev. To be able to operate on digital signal processing systems must have 
a sampling rate of at least 2 x frequency data in accordance with the terms Nyquist, in this design used 200 Hz (fs). Bandpass 
digital filter required specification, the system as follows: 
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• Range band pass filter of 2 Hz to 3 Hz with attenuation that is still allowed no more than 3 dB at a frequency of two 
sides of the ramp. 

• Attenuation should be less than 18 dB at the frequency of 1 Hz and 4 Hz. 


• The type of response required by the ripple Chebyshev 3 dB. 


• Looking folding frequency (fO) is fO = fs / 2 = 200/2 = 100 Hz. There are two frequency of 2 Hz and 3 Hz to be selected 
as a band-edge frequency as fl and f3. For search, the values of the main and sub prime rates, the frequency will be 


changed to a form of normalization, 

Vl = fi/fo = 2/100 =0.02 (9) 

v 3 = f 3 /f 0 = 3/100 =0.03 (10) 

b = v 3 - V! = 0.03 - 0.02 = 0.01 (11) 

v a = f rl /f 0 = 1/100 = 0.01 (12) 

v b = f r2 /f 0 = 4/100 = 0.04 (13) 


The center frequencies v2are: 


7T 7T 7T 7T 71 

tan 2- v 2 =tan — Vi tan — v 3 = tan — x 0.02tan — x 0.03 
2 2 2 2 2 


= tan (l,8)tan (2,7) 
= 0.0314x0.0472 
= 0.0015 

7T 

tan v 2 = 0.0385 

2 


IL v 2 = 2,2048 
2 

v 2 = 0.0245(14) 

Transformation for the constants D and E with the assumption that A* = 1 rad/s. 

D = At cot — = At cot — = cot0.9 = 63.6537 

2/0 2/o 

So to Aa and lb representing two analog frequencies, are as follows: 

A a cos nv 2 — cos nv a 


D 


Aa = D 


sin nv„ 


cos nv 2 — cos nv a 
sin nv n 


Aa = 63.6537 


cos ttO. 0245 — cos ttO.01 


Aa = 63.6537 


sin jtO.01 

cos (4.41) — cos (1.8) 


Aa = 63.6537 


sin (1.8) 

(0.997 - 0.999) 
0.0314 
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Aa = -4.0543 

cos nv 2 — cos nv b 


(17) 


Xb = D 


Ab = 63.6537 


sin nv b 

cos ttO. 0245 — cos7t0.04 


Ab = 63.6537 


sin ttO.04 

cos (4.41) — cos (7.2) 


Ab = 63.6537 


sin (7.2) 
(0.997- 0.992) 


0.1253 

Ab = 2.54 


(18) 


The transfer function of the analog normalized second-order Chebyshev filter with a3d Bripple of Table 2, with the use of 
dummy variables p, that function becomes as below: 


G(p) = 


0.708 


Transformation needed 


0.708 + 0.645 p + p 2 

p = D f(l- Ez -1 + z -2 ) 


P 


(1-z- 2 ) 


p = 63.6537 


r(l- 0.92 z~ 1 + z“ 2 )l 

l (1-z- 2 ) J 


Substituting (20) into the equation(19), 


H(z) = 


G(p) 


p=D 


(1- E Z 1 + Z Z ) 


H(z) = 


0.708 


G(p) 


p =63.6537 


(l - z“ 2 ) J 

(1- 0.92 Z _1 + Z -2 ) 


(1- Z _z ) 


H(z) = 


H(z) = 


0.703+0.645 p+p 2 p =63.6537 

0.708 (1- z~ 2 )z~ 2 


(1- 0.92 Z 1 + Z 2 )1 


(1- Z~ 2 ) 


( 21 ) 


0.708 (1- z~ 2 )z~ 2 + 41,054 (1- z~ 2 ) + 4051,3225 (l- 0.92 z _1 + z~ 2 ) 2 


( 22 ) 
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H(z) = 


0.00017(1- z~ 2 )z~ 2 

1 - 1,8123 z -1 + 2,7894 z" 2 - 1,7940 ;r 3 + 0,96027 ;r 4 


(23) 


the coefficient of the equation (23) are as follows: 


Table 2 

Value Coefficient of Equation (23) 


No 

Koefisien 

Nilai 

Desimal 

Hexadesimal 

1 

A0 

0.00017 

17 

0011 

2 

A2 

-0.00034 

-34 

FFDE 

3 

A4 

0.00017 

17 

0011 

4 

B1 

-1,8123 

-181230 

D3C12 

5 

B2 

2,7894 

278940 

4419C 

6 

B3 

-1,7940 

-179400 

D4338 

7 

B4 

0.96027 

96027 

1771B 



B4 


Figure 7. Realization for Systems of Equations23 
2.7.2 HR filter Deployment Process 

Analog input signal must first sample and digitally using ADC (analog -to-digital con-verter). The resulting binary numbers, 
representing the values of consecutive samples of the input signal, will be transferred to the processors, which perform 
numerical calculations on them. These calculations usually involving mathematical operations by downloading galikan input 
values by a constant and Menam-even the same product. If necessary, result-Tungan perhi, who now represents the values of 
the samples of the filtered signal, is output through a DAC (digital to analog converter) to convert the signal back to its 
original analog form. 

Modules are used for implements its-invest the HR filter is as follows: Module Multiplication, Addition module, splitter 
modules, and modules Reduction, that all of it was developed using the programming language VHDL. 
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2.7.3 Digital Filter Design with MatLab Simulink 

Band-Pass Filter: The specification of the filter as below: 

• Response Type: Band-Pass 

• Design Method: HR Chebyshev 

• Filter Order: 4 th 

• Frequency specification: Fs (200 Hz), Fstopl (2 Hz), Fstop2 (3 Hz) 

• Magnitude Specification: 60 dB (Astop) 

The test results with MatLab program with the input of order 2 are required. Order test follows the tables on MatLab program 
with input from 2, 4, 6, ..., 20. 



Table3. 

Simulation Results Band-Pass Filter with Simulink 


No 

X(n) 

Y(n) Simulink 

1 

2 

0,00200 

2 

4 

0,00400 

3 

6 

0,00600 

4 

8 

0,00801 

5 

10 

0,01001 

6 

12 

0,01202 

7 

14 

0,01402 

8 

16 

0,01602 

9 

18 

0,01803 

10 

20 

0,02003 


The simulation results using MatLab with multiple inputs of different frequencies, among others: 



Figure 9.1nput Frequency 2.3-2.7Hz 
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2.7.4 Implementation of Simulation-Based FPGA Digital Filter 


Band-Pass Filter: coefficient that has been designed (calculated) has been given in the equation above. Now this VHDL code 
used to rate circuit with the use of tool SINTESI Xlink and simulated using the simulator Xlink with the input sequence (2, 4, 
6, 20). At this stage the researchers used a type ofprogram Xlink with Spartan-3E, according to the FPGA board that had 

been prepared. 



SffltTineMliSlOjOOOfS 


Figure II.Simulation Results in FPGA 


2.7.5 Analysis and Discussion 

At the time of testing or simulation is no occurrence of errors or differences in results between MatLab programs with Xilinx. 
The incidence of errors is limited to the types of data types Xilink, about fractions, as in Table4. 


Table 4. 

Simulation Results Band-Pass Filter Based on Simulink with Xilink 


No 

X(n) 

Y(n) Simulink 

Y 2 (n)Xilink 

1 

2 

0,00200 

0 

2 

4 

0,00400 

0 

3 

6 

0,00600 

0 

4 

8 

0,00801 

0 

5 

10 

0,01001 

-1 

6 

12 

0,01202 

-1 

7 

14 

0,01402 

-1 

8 

16 

0,01602 

-1 

9 

18 

0,01803 

-1 

10 

20 

0,02003 

-1 


Band-Pass circuit design has been completed using VHDL and implementation da-lam Xilinx Spartan-3E (package: FG320, 
speed: -5), in the process of making this design using theXilinxISE14.6 design tools. The tables use the following resources 
for Spartan-3E on the design and implementation of digital filters lawyer-data rule, 
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Table 6. 

Advanced HDL Synthesis Report 


Circuit 

Total 

multipli 

ers 

Total 

Addres 

s 

Total 

Registers 

(Flip- 

Flop) 

Low 

Pass 

4 

5 

64 

Band 

Pass 

6 

7 

256 

High 

Pass 

4 

5 

64 


III. Discussion 

The minimum period for a bandpass filter is 1.384 ns (maximum frequency: 722 648 MHz). Number 4 input LUTs for 
bandpass filter is 209. The total number of gates use in the design of implementation is 2408 to bandpass filter. 

IV. Conclusion 

Designs of Digital HR Bandpass Filters have been achieved. A number of findings have 
identified during the design and simulation phases, as follows: 

1 . HR filter can be implemented based FPGAs. 

2. The coefficient of the HR filter design has been done by calculation. VHDL has been used to get into the hardware 
description. VHDL code has been written, synthesized, mapped and then successfully configured with a prototype. 
HR filter is designed completely in accordance with the design requirements. 

3. In testing the filter circuit is still a difference between the results of Matlab simulation program to program Xilinx. 
This is influenced by the use of dividers components that perform the division process with the result of rounding. 

4. One of the problems in the process simulation program Xilinx is the use of integer data, where data from the fraction 
should rounding first. 

5. The results of the FPGA board to produce the same output analog signals with analog signal input. So that this 
algorithm is correct and appropriate, the file * .ucf can be used to make chips infrasound. 
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